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Askaryan Radio Array (ARA), a large-scale radio Cherenkov observatory which scientists propose 
to develop in Antarctica, aims at discovering the origin and evolution of the cosmic accelerators 
that produce the highest energy cosmic rays by means of observing the ultra high energy (UHE) 
cosmogenic neutrinos. To optimize ARA's angular resolution of the incoming UHE neutrinos, which 
is essential for pointing pack to its source, the relation between the reconstruction capabilities of 
ARA and its design is studied. It is found that with the noise effect taken into account, in order 
to make this neutrino angular resolution as good as possible and detection efficiency as high as 
possible, the optimal choice for ARA geometry would be the station spacing of 1.6 km and the 
antenna spacing of 40 m. 



INTRODUCTION 



A limit on the cosmic ray energy was suggested in 1966 
by Kenneth Greisen (US) [1 and Vadim Kuzmin and 
Georgiy Zatsepin (Russia) [2j independently based on in- 
teractions between the cosmic ray and the photons of the 
cosmic microwave background radiation. They predicted 
that cosmic rays with energies over the threshold energy 
of 6x 10 19 eV would interact with cosmic microwave back- 
ground photons to produce pions. This interaction would 
continue until their energies fall below the pion produc- 
tion threshold This theoretical upper limit on the energy 
of cosmic rays from distant sources will create a cutoff 
in the cosmic ray spectrum right at the energy level of 
6 x 10 19 eV. And thus we call this GZK limit or GZK cut- 
off. Furthermore, the interaction of photons and protons 
does not stop at pion productions. These pions continue 
to decay into neutrinos. The whole interaction is named 
as GZK process, and the neutrinos produced from GZK 
process are called GZK neutrinos. The energy of GZK 
neutrino is also very high, still above the level of 10 19 
eV, so it can deserve the name of UHE neutrino. Fig- 
ure [T] shows the energy spectrum of UHECR observation 
and GZK neutrino prediction. The GZK neutrino models 
in this figure were propsed by Kalashev [3], Protheroe, 
and Johnson [4 et al., and UHECR observation data 
are taken from Auger [5], Yakutsk [6 , the Fly's Eye [7], 
AG ASA [8], HiRes [9], and Haverah Park [10]. Error bars 
here only include statistical errors. 
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Figure 1: Energy spectrum of UHECR observation and 
GZK neutrino prediction. 

Since UHE neutrinos can be a proof of the GZK pro- 
cess, but also a key to unveil the mystery of the cosmic 
accelerator, detection of them is scientifically important. 
Neutrinos cannot be directly observed, they but can be 
indirectly observed through their interactions with or- 
dinary matter. There are three possible neutrino inter- 
actions: the elastic scattering, the charge current (CC) 
interaction, and the neutral current (NC) interaction. 
The outgoing high energy particles would, due to the 
Askaryan effect [IT], result in 20% excess of fast moving 
negative charges at the shower maximum and produce 
Cherenkov radiation, the radio band of which is coherent 
in ice and can be employed as a probe for UHE neu- 
trinos. Askaryan effect states that a high energy parti- 
cle which travels faster than the light in dense dielectric 
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material can lead to charge asymmetry in its shower be- 
cause positrons in the shower have higher probability of 
being stopped by the atoms in the medium and in ad- 
dition the bound-state electrons tend to be ionized and 
comove with the shower. Then these fast moving ac- 
cess charges would cause Cherenkov radiation, where the 
radio and microwave bands of which are coherent. By 
invoking the Askaryan effect as a mechnism for detec- 
tion, UHE neutrinos could in principle be observed. One 
of such attempts is the Antarctic Impulse Transient An- 
tenna (ANITA) project, which intends to detect cosmo- 
gentic neutrinos interacting with the Antarctic ice. The 
ANITA Collaboration has performed an experiment at 
the Stanford Linear Accelerator Center (SLAC) in June 
2006, which confirmed this effect in ice [12]. 

Askaryan Radio Array (ARA), a large-scale radio 
Cherenkov observatory which scientists propose to de- 
velop in Antarctica, aims at discovering the origin and 
evolution of the cosmic accelerators that produce the 
highest energy cosmic rays by means of observing the 
ultra high energy (UHE) cosmogenic neutrinos. The rea- 
sons why ARA choose the Antarctic as the experiment 
site are the following: 1. There is plenty of ice as the 
target for detecting neutrinos. 2. The ice is so trans- 
parent to the RF shower signal that the spacing of RF 
detectors can be sufficiently large to enhance the effec- 
tive volume and event rate. 3. It is more radio-quiet 
than other places in the world so as to reduce artificial 
signals considerably. 4. The temperature is so low that 
the background noise also reduces considerably. 

For these reasons, Antarctica is a very proper site to 
do the UHE neutrino experiment. With such a nature- 
given experimental environment, the next issue would be 
how to optimize the array geometry so as to maximize 
the performance. The primary goal of this paper work 
is to assess and optimize the capability of ARA, particu- 
larly the capability of reconstructing neutrino incoming 
directions, by means of Monte Carlo simulations. 



II. SIMULATION METHOD 
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Figure 2: The geometry of the distribution and their 
coordinates. 
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Figure 3: ARA antenna cluster geometry in a station, 
where there are twelve antennas, indicated by green 
squares. 



A. Setting Array Geometry 

The proposed radio-based neutrino detector array, 
ARA, will eventually cover about 80 km 2 at the South 
Pole. There will be 37 antenna stations in the complete 
ARA. These 37 stations are located on a hexagonal lat- 
tice, as shown in Fig. [2j with a station spacing of 1.33 
km. Note that the coordinate in this figure and in this 
analysis has its origin defined at the center of ARA, on 
the surface of the ice, and the z axis points to the sky. 

Each station is composed of a triad of boreholes with 
depths of 200 m, on the corners of an equilateral tri- 
angle. Each borehole has four antennas, two of which 
are the horizontal-polarization (Hpol) antennas and the 
other two the vertical-polarization (Vpol) antennas, as 
shown in Fig. [3j A pair of antennas, a Hpol one and 



a Vpol one, can detect the strengths of electrical field 
projected to horizontal plane (2D) and vertical line (ID) 
respectively, and the find the possible direction of the 
electric field. The side length of the equilateral triangle 
and the distance between two Vpol antennas in a bore- 
hole, are set the same, at 30 m. The location coordinate 
of antenna i is denoted as x.f nt . 



B. Event Production 

In the simulation, the shower events resulted from the 
CC or NC interactions are generated in the ice. This 
event generation does not differentiate neutrinos or anti- 
neutrinos, and flavors are not considered, either. In 
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each event, 6 parameters are generated, including the 
shower location, #f ^ n ? 2/f ^ n ? the moving direction 

of neutrino: gen , (\) g v en , and the intensity of the radio 
Cherenkov radiation induced by the shower followed by 
this interaction, V^ en . We treat the shower location the 
same as the neutrino interaction vertex because of the 
small shower size in ice. The generated shower locations 
are uniformly distributed over a cylinder volume, where 
the center of the cylinder volume is located at the cen- 
ter of ARA. This volume has an axis along the vertical 
direction passing through the center of ARA, and has a 
radius of 6 km and a height from z = -2 km to z = 0. The 
reason why we set the height of this event cylinder as 2 
km is that the thickness of the ice in Antarctic is approx- 
imately 2 km. The choice of 6 km for the radius is due 
to the following reason. For the events with shower lo- 
cations far away from ARA detectors, its radio radiation 
can not reach ARA because of attenuation. Therefore, 
the farthest distance for the radio signal of the event to 
travel to the ARA center is approximately estimated as 
lkm + 1.33km x 3 = 5km. For safety reason, we set it 
as 6 km rather than 5 km. Black dots are the generated 
shower locations, distributed uniformly in this field. Red 
circles represent the 37 stations. 

The moving directions of the neutrinos are also set 
uniformly distributed isotropically over 4 tt solid angle. 
Furthermore, the outcome of a recorded waveform has 
been converted into voltage from electric field through 
the readout electronics and thus we set the initial inten- 
sity of the Cherenkov radiation in terms of the voltage, 
V£ en , in the range of to 5 V. 

At this stage, we generate six parameters for each 
event: *^ n , 0g ew , V gen . The first three 

are also denoted as the shower location vector, x^ n . The 
next two can also be described by a unit vector, p^ en . 
There are 300 simulation events generated in each simu- 
lation. 

And these six event parameters are to be determined 
through x 2 fit in the next few sections. 



C. Radio Cherenkov Wave from the shower 
location to Antennas 

The Cherenkov radiation is set as a point source radia- 
tion because the shower size is of the order of ~m and the 
propagation length before being received is ~km. The 
radiation wave front has a cone shape with the apex at 
the interaction location, with the axis along the neutrino 
moving direction, and the span angle of the RF wave 
from 55 degree to 57 degree. The voltage waveform of 
this radiation signal is set as a bipolar wave: 



y = x . e -* 2 /^\ 



where y is voltage and x stands for time. 



(i) 



As this wave propagates through the ice, the wave in- 
tensity changes as 
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cosaf en (for Vpol antenna), 



where .Do is the distance from the shower location to 
the location where V§ en is measured, 1 km, 6? en is the 
separation angle between the vector pg en and the vector 
x ant _ x ^ n , and af en is the separation angle between the 
direction vector of the antenna z, (x, y, z) = (0, 0, 1) and 
the direction vector of the electric field, (x^ n£ — x^ n ) x 
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The travel time of this signal is 
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where c is the speed of light divided by the refraction 
index of ice. 

At the signal receiving end, the oscilloscope has time 
bin of 0.39 ns, and the time window is 100 ns. Noise be- 
fore circuit has Gaussian distribution with mean voltage 
of and cr n oise 0.035 mV, whereas the trigger thresholds 
are that the Cherenkov cone intersects with the antenna 
and the attenuated signal must be larger than la noise- 
An waveform is made by the following steps: an orig- 
inally bipolar waveform magnified by a factor of signal 
strength V[ ea \ shifted to the right by a time lag of ^, 
and then added with noise. 



D. Determination of Arrival Time Difference and 
Pulse Voltage 

To do the reconstruction of the events in the next stage, 
we have to extract arrival time difference, At° 6s , and the 
pulse voltage, V° bs , from the waveform of each triggered 
antenna. 

The time when the signal arrives at the antennas 
should be precisely determined, and then with the dif- 
ference of arrival time between any two antennas, and 
with the requirement that at least four antennas must be 
triggered, the shower location can be obtained through 
the process of fitting. One way can be applied to cal- 
culate arrival time, t? 6s , for each antenna is the use of 
the point where V = between the maximum and the 
minimum amplitudes. 

From the procedure described in the previous para- 
graph, for each antenna we can obtain an arrival time. 
Arrival time difference, which is the information actually 
used in the reconstruction, is the arrival time subtracted 
by the reference of the arrival time, £q 6s , which is defined 
as the arrival time for the antenna receiving the strongest 
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signal among all antennas. Therefore, 
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(4) 



As for the pulse voltage, it is either the maximum point 
or the minimum point, depending on which one arrived 
at the antenna first. 



E. Reconstruction of Neutrino Moving Directions 

Our event reconstruction procedure is divided into two 
stages. The first stage is the reconstruction of shower 
location, x^ n . In this stage, the needed information is 
the arrival time difference, At° bs for each antenna. We 
set up a x 2 formula: 
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Figure 4: Resolution of shower location in z axis. 
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(5) 



where i is the index for all the triggered antenna, and 
At^ yp is the hypothesized arrival time difference. By 
minimizing xti the best fit x^ c can be found, where a 
grid search is employed. Local minima of xl value in 
the hypothesized- variable space is a serious problem and 
prohibit us from using other efficient ways to find the 
global minimum. 

In the second stage of reconstruction, we still use x 2 
to find the best-fit. In this stage, the moving direction of 
neutrinos, p^ en , is to be reconstructed, and the needed 
information is the pulse voltage received in each antenna. 
Furthermore, we also have to input the reconstructed 
shower location, x^ c , which is obtained in the first stage 
of reconstruction. Otherwise, we have to treat it as an 
unknown parameter to be reconstructed and this would 
intensively increase the computing time. The x 2 formula 
in the second stage is given as 
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obs 



vl iyp (x r s e h c ,Vo vp .u h yp)V 



(6) 
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where i is the index for all the triggered antenna, and 
V- hyp is the hypothesized pulse voltage. By minimizing 
(l,6£ ec ,6£ ec ) can be found. 



xl? the best fit 



III. RESULTS 

A. Resolutions of Shower Location, RF Wave 
Direction, and Neutrino Moving Direction 

With Eq. [5j the shower location can be reconstructed. 
The obtained resolutions of the shower location are 0.143 
km in x axis, 0.098 km in y axis, 0.07 km in z axis. 
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Figure 5: Resolution of RF wave direction in zenith 
angle. 



The resolution is the RMS value in Fig. |4j which is the 
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distributions of Az 

Once the reconstructed shower locations are obtained, 
these reconstructed locations are taken as input in Eq. [6] 
for the reconstruction of neutrino moving directions. 

The obtained resolutions of the reconstructed RF wave 
direction, are 1.45° in direction , and 3.69° in direc- 
tion. The resolution in direction is shown in Figs. [5j 
which is the distributions of AOrf = Q'rf — @rf • ^ ne ^ 
and (j) here are the zenith angle and the azimuthal angle 
of the spherical coordinate with origin defined as the lo- 
cation of the antenna receiving the strongest signal, and 
z axis as before. 

After the shower location and the RF wave direction 
are obtained, one can compute the neutrino direction ac- 
cording to Eq. [6] The obtained resolutions, in 9 direc- 
tion 4.88°, and in (j) direction 3.76°, which are the RMS 
values of the distributions of A9 V = 9 r v ec — #g en , and 
A<p v = (jf v ec — <\) 9 v en . The average of the separation angle 
between the generated v direction and the reconstructed 
one is shown in Fig. [6| which is 2.38°. 

The average of this angle is taken for the comparison of 
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Figure 6: Distribution of the separation angle between 
the generated v direction and the reconstructed one. 
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Figure 7: Resolutions of neutrino direction, (AQ„) 
versus antenna spacings and station spacings. 



the neutrino angular resolution in this analysis because 
the separation angles are always positive, and thus the 
RMS value may not represent a proper indication of res- 
olution. Note that and (j) here are the zenith angle and 
the azimuthal angle of the spherical coordinate with the 
origin defined as the location of ARA center, and z axis 
as before. As mentioned before, the results presented so 
far have employed the ARA array geometry of station 
spacing as 1.33 km and antenna spacing as 30 m. 
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B. Optimization of ARA 

To optimize the ARA, 16 different antenna spacings 
and 10 different station spacings are selected for the study 
on the resolution of the neutrino moving direction and the 
detection efficiency along with studies of noise effect. The 
optimum would be achieved when the resolution of the 
neutrino moving direction, i.e. (AO^), is as good as pos- 
sible, and the detection efficiency is as high as possible. 
The detection efficiency is defined as the number of trig- 
gered events that pass the trigger threshold divided by 
the total number of generated events in the cylinder vol- 
ume, where the threshold applied to the pulse voltage is 

7 O noise- 

The antenna spacing varies from 10° 7 m to 10 2,2 m 
in steps of 0.1 in the power index of 10. The station 
spacing changes from 1.33km/5 to 1.33km x 2 in steps 
of 1.33km/5. Note that the antenna spacing means the 
distance from the top antenna to the bottom one. The 
vertical spacings between any two antennas are the same, 
and the center of the four antennas in a borehole is lo- 
cated at the depth of 200 m. In addition, the side of 
the equilateral triangle in a station is set the same as the 
antenna spacing. 

The mean value of the separation angles (AQ U ) versus 
the antenna spacings is shown in Fig. [7| whereas the de- 
tection efficiencies versus the antenna spacings are given 
in Fig. H 



Figure 8: Detection efficiencies versus antenna spacings 
and station spacings. 



Figs. [7| and [8] suggest that (AQ U ) can be less than 5° 
if the station spacing is set in the range of 1.33 km to 1.9 
km and the antenna spacing is set in the range of 40 m 
to 100 m. One may notice that the detection efficiency 
reach a saturated value, ~ 70%, when the station spacing 
is grater than ~ 1.5km. 

To finalize the optimal choice for the ARA geometry, 
the effects of different noise levels added to the origi- 
nal waveform and different trigger thresholds are stud- 
ied, too. The value of cr no i se is set at 0.035 mV for all 
analysis presented so far with V^ en varying in the range 
of to 5 V. In the following studies of how the noise 
levels would affect the resolution of the neutrino mov- 
ing direction, in each case a different level of noise added 
to the waveform is assumed, i.e. cr noise = 06a no i se , with 
a greater than one, whereas V^ en is fixed at 5 V. Dif- 
ferent trigger thresholds are applied: V° bs > 3.5a noise , 
V° bs > 7a no i se . For these studies, only 100 events are 
generated in each case. The results of (A6„) and the 
detection efficiency versus the noise level under different 
trigger thresholds are presented in Figs. ?? to ?? for dif- 
ferent antenna spacings and different station spacings. It 
was found that the larger cr' noise added to the waveforms, 
the worse the resolution of the neutrino moving direction, 
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Figure 9: Resolutions of neutrino direction, (AG^), 
versus different noise levels, and antenna spacings, 
where station spacing is set at 1.60 km and the trigger 
threshold is 3.5cr no i se . 



Figure 11: Detection Efficiencies versus different noise 
levels and antenna spacings, where station spacing is set 
at 1.60 km and the trigger threshold is 3.bcr no i se . 
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Figure 10: Resolutions of neutrino direction, (AG^), 
versus different noise levels, and antenna spacings, 
where station spacing is set at 1.60 km and the trigger 
threshold is lcr noise . 

which is as expected. In addition, the higher the trigger 
threshold, the lower the detection efficiency. 

In summary, with the noise effect taken into account, 
in order to make the resolution of the neutrino moving 
direction as good as possible and detection efficiency as 
high as possible, the optimal choice for ARA geometry 
would be 1.6 km for the station spacing and 40 m for the 
antenna spacing. 



IV. SUMMARY 

Angular Resolution of Neutrino Moving Direction: 
One of the main goals of ARA is to point back to cos- 
mic accelerators through the determination of the UHE 
neutrino moving directions, so the resolution of it is par- 
ticularly important. 



Figure 12: Detection Efficiencies versus different noise 
levels and antenna spacings, where station spacing is set 
at 1.60 km and the trigger threshold is 7cr noise . 



To optimize the ARA, both the resolution of the neu- 
trino moving direction and the detection efficiency should 
be considered. Basically, the detection efficiency in- 
creases as the station spacing gets larger. From Fig. [§J 
however, it reaches a plateau of ~ 70% detection effi- 
ciency when the station spacing is grater than ~1.5 km 
where the regions which each station can cover no longer 
overlap. With the noise effect taken into account, in order 
to make the resolution of the neutrino moving direction 
as good as possible and detection efficiency as high as 
possible, the optimal choice for ARA geometry would be 
1.6 km for the station spacing and 40 m for the antenna 
spacing. 

In the simulation of angular resolution of neutrino di- 
rection for Antarctic Ross Ice Shelf ANtenna Neutrino 
Array (ARIANNA) experiment, the resolution in di- 
rection is 1.1° [15]. However, to reach such a good reso- 
lution, ARIANNA has to build its array up to 11 stations 
per km 2 , which means that its antenna density has to be 
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13 times greater than ARA if we set the station spacing 
as 1.33 km. Based on this comparing, the design of ARA 
is in a better balance point between the resolution and 
the cost. 

In the future, if ARA can get more funding to increase 
the density of the antenna number, a much better reso- 
lution of neutrino moving direction can be achieved. 
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